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Introduction 

The National Park Service Protective Glazing Study conducted by Inspired Partnerships 
has provided a unique opportunity for the quantitative analysis of the effect of 
protective glazing. Dataloggers were installed in windows with protective glazing at St. 
John U. C. C. in Evanston. Only a small portion of the available data has been 
analyzed. Preliminary analysis has concentrated on data from two similar windows with 
the same geometry, both with protective glazing, but one with venting. This gives a 
direct measure of the effect of venting. 

Analysis 

The first requirement of this analysis was the calculation of the solar input (insolation) 
for any location on the earth, any date, any time, and any orientation (elevation and 
azimuth) of the surface receiving the solar input. Exact calculations would require 
orbital mechanics, but a good approximation was found which was programmed in the 
HP keystroke language and had an associated writeup of the equations. These 
equations were then programmed in Visual Basic as an Excel function. It is listed as 
"radvmod" in the appendix. The first check was to compare calculations of the Excel 
function with calculations with the keystroke program on an HP calculator. After 
correcting several errors in the writeup of the equations, satisfactory agreement was 
obtained. The second check was to compare the Excel function, listed in the appendix 
as "radv", output with tables of solar insolation accepted in the field. Comparisons were 
made with data in tables from "Solar Heating and Cooling", Jan F. Krieder and Frank 
Kreith (Ref. 1). Figures 1 through 3 present comparisons on March 21 at 40 deg N Lat, 
for vertical, 40 deg inclination, and horizontal surfaces. Very good agreement is 
shown. 

Figures 4 through 6 the same comparisons for June 21 . The vertical surface 
comparison, Figure 4, is quite interesting. An E-W wall in summer does not receive 
direct sunlight for some time after sunrise or for some time before sunset. The tables in 
Ref 1 include ground reflection, so some insolation occurs before the wall gets direct 
sun. This is shown graphically in Figure 4, since the approximate calculation does not 
include ground reflection. Figures 5 and 6 show that the calculation is a little high in 
midday, but is still quite acceptable. 

The calculations for September 21 are the same as those for March 21, so the 
comparisons are not repeated for September. 

Figures 7 through 9 present results for December 21 , the most southward travel of the 
sun. Calculations are just slightly below the table values for a vertical wall. 
Calculations for the inclined and horizontal surfaces are not quite as good, as much as 
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14 per cent low at noon, but are still acceptable. Fortunately, the primary interest in this 
application is in the vertical surface. 

The second requirement of this analysis is a solution for the coupled differential 
equations of the heat transfer in the window and cover. Some of the solar input is 
reflected from the cover, some is adsorbed by the cover, another portion is absorbed 
by the lead came and the glass of the window, with the remainder adsorbed by the 
interior structure. 

The first step in satisfying the second requirement is to evaluate the physical processes 
in as complete a manner as possible. First, some of the solar input is reflected from the 
cover. Reflection from a transparent medium is a function of the refractive index and 
the angle of incidence of the solar input. The Handbook of Chemistry and Physics 
(Ref. 2) gives a table on p. 2437 for reflection for a refractive index of 1.55 (appropriate 
for both glass and Lexan). Figure 10 shows that a simple exponential equation is a 
good representation of the results calculated using Fresnel's formula. Note that the 
solar input has a constantly changing angle with respect to the cover, so an equation is 
required. 

Second, an estimate of the solar energy adsorbed by the Lexan is required. A phone 
call to GE Plastics Structured Products Technical Sales produced a fax of a test report 
on the exposure of UV-coated and uncoated Lexan to the sun (Ref. 3). Figure 1 1 
presents the data from the fax and an equation which is a reasonable fit to the data. 
An extrapolation of the data on uncoated Lexan to 14 years (exposure time of these 
covers), accounting for reflection, gave a value of 0.20 for adsorption of the solar input 
by the cover. A rough check of transmission with a light meter was consistent with this 
value, setting one of the constants in the equations. 

Third, an estimate of the solar energy adsorbed by the stained glass is required. First, 
from the average length in one square foot of the window, about 20% of the area is 
covered by lead, adsorbing essentially all of the solar input. A rough check of 
transmission through several samples of medium density stained glass gave an 
average of about 25% transmission. The effect of reflection and the complete 
adsorption of the lead are approximately in balance, so 75% was taken for the 
adsorption of the stained glass. Note that light meter measurements on individual 
windows would be desirable. Two thermal inertia constants were computed, based of 
the mass and specific heat of one square foot of the Lexan cover (0.55 BTU/deg F) 
and one square foot of the window ( 0.36 BTU/deg F). Other constants required were 
determined by the geometry of the windows^(thickness of air space, window area, vent 
area in and out, vertical distance between the vent areas), the properties of air, and the 
gravitational constant. Average atmospheric pressure was given by weather reports for 
O'Hare Airport (Ref. 4). 

With the best estimates or determination of the constants complete, the next step is to 
select the heat transfer equations to be applied to the components of the window and 
cover system. First, the heat transfer from the interior to the stained glass window is 
driven by natural convection. The flow velocity across the window is a function of the 
temperature difference between the window and the interior. Marks "Mechanical 
Engineers' Handbook (Ref. 5, p 374) gives an expression for natural convection for a 
vertical surface over one foot high (Eq. 1). 
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U a = natural convection heat transfer ~ BTU I ft 2 / deg F 
T f = inside temperature ~ deg F 
T; = stained glass window temperature ~ deg F 
[7^-7^|= absolute value of the temperature difference ~ deg F 

Then the heat transfer from the room to one square foot of the window is given by 
Equation 2. 

Q>=&*<T t -T t ) (2) 

Q u - heat transfer ~ BTU I ft 2 Ihr 

Note that the equation is written in terms of heat transfer to the window. That is, if T t 
the interior temperature, is hotter than T s , the window temperature, the window is 
gaining heat from the interior. If the window temperature is hotter than the interior 
temperature, the window is losing heat to the interior. 

The next step is the heat transfer from the window to the air space to the cover. Since 
the thermal inertia of the air in one square foot of air space is very small, the heat flow 
across the air space can be assumed to be quasi-steady. That is, the temperature of 
the air space is approximately equal to the average temperature of the window and the 
cover. Small transient deviations may be produced by venting. "Heat Loss Through 
Stained Glass Windows", by Wayne Simon (Ref. 6) presents the variation in resistance 
factor for enclosed air spaces for depths up to three quarters of an inch, and states that 
the resistance is constant for larger depths. As a convenience for this program, this 
data has been fitted with an exponential equation, Figure 12, so that the program 
constant can be just the depth of the air space. Now the heat transfer from the window 
to the cover is, except for the effect of venting, equal to the heat transfer from the 
window to the air space and to the heat transfer from the air space to the cover. From 
Figure 12, Equation 3, for R x , the resistance to heat flow from the stained glass 
window to the cover, can be written. 

R x = 0.83(1 -e 6iD -) (3) 

D a =air space depth -in 

Then, since the U-factor is just the reciprocal of the R-factor, 

And 

Q 1C = (T S -T C )-U SC (5) 
Q x = heat from s to c~ BTU I ft 2 Ihr 



Therefor, the heat transfer from the air space to the stained glass window can be 
written: 

Q as =2-{T a -T s )U sc (6) 

The last input to the stained glass window is the portion of the solar radiation adsorbed 
by the window, this is given in Equation 7. 

RAD S = {solar insolation) ■ (1 - reft) ■ (1 - cover adsorb) ■ {window adsorb) (7) . 
Then the differential equation for the temperature of the window can be written: 

dt TM S 

TM S = thermal mass of window ~ BTU I ft 2 1 deg F 

The next task is to develop the differential equation for the temperature of the air 
space. The first two terms are simple, since the heat transfer to the air space is just the 
negative of the heat transfer from the window and the cover to the air space, as shown 
in Equations 9 and 10. 

Q sa =-Q as = 2{T s -T a )U sc (9) 

Q ca =-Q ae -2-{T c -T a )-U sc (10) 

The nest step is to develop the equation for the effect of venting. When the air in the 
air space is hotter (or colder) than the outside air, the difference in density creates a 
buoyancy force which, if venting area is provided at the top and the bottom of the air 
space, causes a flow of air. Marks Handbook (Ref. 5, p 1 120) gives an equation for the 
draft (pressure difference) as a function of height (vertical distance between inlet and 
outlet vent area), ambient pressure, outside temperature, and the air space 
temperature. Converting to the units of this work, it is given as Equation 1 1 

D = 0.0m-P H \ — — — — 1 I (11) 
^+460 T a +460 J 

D = draft 

P - ambient pressu ~ 

H = vertical distance between inlet and outlet ~ ft 
T a = ambient temperature- deg F 



= air space temperature- deg F 



Note that the quantity "460" added to the temperatures converts to an absolute 
temperature scale, degrees Rankine. Now the steady state flow through the air space 
represents a balance between the draft and the pressure loss due to flow. For the 
usual case, where the area of the vents is much smaller than the flow area through the 
air space, the loss can be evaluated as a fraction of the dynamic pressure through the 
inlet and outlet, as given in Equation (12). 



D = 




n = 


loss ratio 


p = 


density ~ -^r 


G = 


gravitational acceleation- -~ 
secr 


Vj = 


vent inlet velocity ~ — 
sec 


K = 


vent outlet velocity - — 
sec 



Using mass conservation to express vent velocity ratio in terms of vent area ratio, 
Equation 13 can be developed. 



r..£: ' J| ~ (13) 

Note that, in order to compute correctly if the air space temperature is colder than the 
outside air, the sign of D has been moved outside the square root. This could occur in 
an air conditioned building with small solar input. The rate of heat loss (or gain) due to 
the vent flow can now be computed with Equation 14. 

^ jm-p-AVrit-r.) (U) . 

Q vent = venting heat loss ~ BTU I ft 2 Ihr 
= area of stained glass window - ft 2 



Now the differential equation for the temperature of the air space can be written as 
Equation 15. 



dL = Q. + <L+Q~ (15) 

dt TM a 

TM a = thermal mass of air ~ BTU I ft 2 1 deg F 

The last step is to develop the differentia! equation for the temperature of the cover. 
Reference 6 presented a table of R-Factor for the exterior of a building as a function of 
wind speed. For this program, wind speed is a variable, so an equation is required. 
Figure 13 presents a comparison of Equation 16 with the data from Reference 6. 

R oc = 068 (16) 
°° (1 + 0.2-V) 

V - wind speed ~ mph 

Again, since the U-factor is just the reciprocal of the R-factor, 

u -~£ (17) 

And Q ee =(T 0 -T c ) : U oe (18) 

Q oc = heat from outside to cover ~ BTU I ft 2 1 hr 

The heat transfer from the air space to the cover, Q ac , has already been developed in 
Equation 10, so the remaining term, the portion of the solar radiation adsorbed by the 
cover, is given by Equation (19). 

RAD C - {solar insolation) • (1 - re/7) • (cover adsorb) (1 9) 

Now the differential equation for the temperature of the cover can be written as 
Equation 20 

dT c ^ Q ae + Q oc +RAD c . 
dt TM C 

TM C = thermal mass of cover ~ BTU I ft 2 1 deg F 

Note that the three differential equations, (8,15, and 20) are not independent. Each 
temperature appears in the differential equation of the others, making this a coupled 
set of differential equations. The solution is obtained by a very fine time scale 
simultaneous integration (12 second interval) with ten repetitions of the complete 
integration. This Excel function is listed in the appendix as "Tinteg". 



The temperature time histories of the components then depend on the natural 
convection from the window to the interior of the building, the mixed 
conduction/convection between the window and the air space and between the air 
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space and the cover, the forced convection from the cover to the exterior atmosphere, 
and, for the case of venting, the removal of heat from the air space by buoyancy driven 
convection. In actuality, the temperatures also vary with vertical and horizontal location 
in the window and cover, in order to simplify the analysis, it is assumed the 
temperature at all locations in the window and cover is uniform. Another way of stating 
this restriction is that the analysis calculates the space averaged window, airspace, and 
cover temperatures. 

The third requirement for tine analysis is a knowledge of the local weather. First, the 
heat transfer from the cover to the outside atmosphere varies by a factor of four as the 
wind speed varies from 0 to 15 mph, a very modest wind. For this purpose, data for 
O'Hare airport was obtained from the National Climactic Data Center. No direct 
measurement of solar insolation was found, but the weather descriptions were used to 
adjust the variation in transmission through the atmosphere (cloud cover, fog, blowing 
snow, etc.). 

First, the data was searched to find a day as clear as possible. January 26. 1996 had a 
little fog in the morning, but was otherwise clear. This made possible a test of the 
program with minimum ambiguity. January 27 was added as an example of a day with 
much smaller solar input. Figure 14 compares the measured air space temperature of 
two similar windows in the same wall, differing only in venting. One is sealed, the other 
has 0.147 square inch vents at the top and bottom of the cover. The effect of venting 
is not large, but the peak temperatures are significantly reduced. This probably 
indicates that the vent area is on the small side. Larger vent areas will be investigated 
with the computer model. 

Figure 15 makes the same comparison for the measurements of humidity. The vented 
window has significantly lower humidity during solar insolation, even though the 
humidity levels were fairly low in both windows. The equations for humidity have not 
yet been included in the computer model, so the effect of larger vent areas on humidity 
cannot yet be computed. 

Figure 16 compares measured and calculated air space temperatures for the unvented 
window. Figure 17 presents the adjusted air transmission factor which was used to 
compute the calculated temperature. No independent measure of the solar insolation 
was available, so the magnitude of the dips in the transmission factor have been 
determined by the measured air space temperature. The location in time of the dips is 
consistent with the weather reports. Figure 18 presents the measured wind speed at 
O'Hare airport, which is also a vital part of the calculation through its effect on exterior 
heat transfer. 

Comparison of measured and calculated air space temperatures for the vented window 
is given in Figure 19. Note that the adjusted air transmission factor is that of Figure 17 
and, of course, the measured wind speed is that of Figure 18. The agreement is not 
quite as good as that for the unvented window. This may indicate that significant 
vertical variation in the air space temperature may occur with venting. 

Figure 20 presents the computed effect of increasing the vent area. It is seen that the 
vent area of the test was too small to provide significant reduction of temperature from 
that of an unvented window. . 
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Conclusions and Recommendations 



Preliminary analysis of the data is quite encouraging. The only missing factor is an 
independent measurement of solar insolation. The humidity calculations can be added 
to the program at a later time. A real test of the humidity calculations will require data 
from the southeast region, where air conditioned interiors and a hot humid exterior 
atmosphere create real humidity problems. Finally, for a reliable prediction of the effect 
of protective glazing on stained glass windows, the computed results must be 
compared with measured data for several different installations. The installations 
should be chosen with consideration of the locations for which the National Climactic 
Data Center has detailed weather records. 
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Figure 1 - Solar Insolation in BTU/sq ft/hr 
(March 21, 40 deg N Lat, Vertical E-W Wall) 
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Figure 2 - Solar insolation in BTU/sq ft/hr 
(March 21, 40 deg N Lat, E-W wall at 40 deg Inclination) 
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Figure 3 - Solar Insolation in BTU/sq ft/hr 
(March 21. 40 deg N Lat, Horizontal Surface) 
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Figure 4 - Solar Insolation in BTU/sq ft/hr 
(June 21, 40 deg N Lat, Vertical E-W Wall) 
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Figure 5 - Solar Insolation in BTU/sq ft/hr (June 
21, 40 deg N Lat, E-W Wall at 40 deg Inclination) 
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Figure 6 - Solar Insolation in BTU/sq ft/hr 
(June 21, 40 deg N Lat, Horizontal Surface) 
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Figure 7 - Solar Insolation in BTU/sq ft/hr 
(December 21, 40 deg N Lat, Vertical E-W Wall) 
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Figure 8 - Solar Insolation in BTU/sq ft/hr 
(Oecenber 21, 40 deg N Lat, E-W Wall at 40 deg Inclination 




Figure 9 - Solar Insolation in BTU/sq ft/hr 
(December 21, 40 deg N Lat, Horizontal Surface) 
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Figure 12 - R-Factor for Enclosed Air Space 
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Figure 14 - Measured Air Space Temperature 
(unvented and vented) 





